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Abstract. This article is devoted to the analysis of a Monte-Carlo method to 
approximate effective coefficients in stochastic homogenization of discrete ellip- 
tic equations. We consider the case of independent and identically distributed 
coefficients, and adopt the point of view of the random walk in a random 
environment. Given some final time t > 0, a natural approximation of the 
homogenized coefficients is given by the empirical average of the final squared 
positions rescaled by t of n independent random walks in n independent envi- 
ronments. Relying on a new quantitative version of Kipnis-Varadhan's theorem 
(which is of independent interest), we first give a sharp estimate of the error 
between the homogenized coefficients and the expectation of the rescaled final 
position of the random walk in terms of t. We then complete the error analysis 
by quantifying the fluctuations of the empirical average in terms of n and t, 
and prove a large-deviation estimate. Compared to other numerical strategics, 
this Monte-Carlo approach has the advantage to be dimension-independent in 
terms of convergence rate and computational cost. 

Keywords: random walk, random environment, stochastic homogenization, 
effective coefficients, Monte-Carlo method, quantitative estimates. 

2010 Mathematics Subject Classification: 35B27, 60K37, 60H25, 65C05, 
60H35, 60G50. 



1. Introduction 

This article is part of a larger program, which consists in devising and quanti- 
tatively analyzing numerical methods to approximate effective coefficients in sto- 
chastic homogenization of linear elliptic equations. More precisely we tackle here 
the case of a discrete elliptic equation with independent and identically distributed 
coefficients (see however the end of this introduction for more general statistics), 
and present and fully analyze an approximation procedure based on a Monte-Carlo 
method. 

A first possibility to approximate effective coefficients is to directly solve the 
so-called corrector equation. In this approach, a first step towards the derivation of 
error estimates is a quantification of the qualitative results proved by Kiinnemann 
|Kii83j (and inspired by Papanicolaou and Varadhan's treatment of the continuous 
case |PV79p and Kozlov |Ko87j . In the stochastic case, such an equation is posed on 
the whole Z d , and we need to localize it on a bounded domain, say the hypercube Qr 
of side R > 0. As shown in a series of papers by Otto and the first author [GOlOal 
IGOIObj , and the first author [GUP] , there are three contributions to the L 2 -error in 
probability between the true homogenized coefficients and its approximation. The 
dominant error in small dimensions takes the form of a variance: it measures the 
fact that the approximation of the homogenized coefficients by the average of the 
energy density of the corrector on a box Q r fluctuates. This error decays at the rate 
of the central limit theorem Rr d in any dimension (with a logarithmic correction 
for d = 2). The second error is the so-called systematic error: it is due to the 
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fact that we have modified the corrector equation by adding a zero-order term of 
strength T _1 > (as standard in the analysis of the well-posedness of the corrector 
equation). The scaling of this error depends on the dimension and saturates at 
dimension 4. It is of higher order than the random error up to dimension 8. The 
last error is due to the use of boundary conditions on the bounded domain Qr. 
Provided there is a buffer region, this error is exponentially small in the distance 
to the buffer zone measured in units of yT, 

This approach has two main drawbacks. First the numerical method only con- 
verges at the central limit theorem scaling in terms of R up to dimension 8, which 
is somehow disappointing from a conceptual point of view (although this is already 
fine in practice). Second, although the size of the buffer zone is roughly independent 
of the dimension, its cost with respect to the central limit theorem scaling dramat- 
ically increases with the dimension (recall that in dimension d, the CLT scaling is 
R~ d , so that in high dimension, we may consider smaller R for a given precision, 
whereas the use of boundary conditions requires R 3> VT in any dimension). Based 
on ideas of the second author in |Moll] , we have taken advantage of the spectral 
representation of the homogenized coefficients (originally introduced by Papanico- 
laou and Varadhan to prove their qualitative homogenization result) in order to 
devise and analyze new approximation formulas for the homogenized coefficients in 
[GMlOj . In particular, this has allowed us to get rid of the restriction on dimen- 
sion, and exhibit refinements of the numerical method of [GllOj which converge at 
the central limit theorem scaling in any dimension (thus avoiding the first men- 
tioned drawback). Unfortunately, the second drawback is inherent to the type of 
method used: if the corrector equation has to be solved on a bounded domain Qr, 
boundary conditions need to be imposed on the boundary OQr. Since their values 
are actually also part of the problem, a buffer zone seems mandatory — with the 
notable exception of the periodization method, whose analysis is yet still unclear to 
us, especially when spatial correlations are introduced in the coefficients. In order 
to avoid the issue of boundary conditions, we adopt here another point of view on 
the problem: the random walk in random environment approach. This other point 
of view on the same homogenization problem has been analyzed in the celebrated 
paper [KV86] by Kipnis and Varadhan, and then extended by De Masi, Ferrari, 
Goldstein, and Wick [DFGW89 . The strategy of the present paper is to obtain 
an approximation of the homogenized coefficients by the numerical simulation of 
this random walk up to some large time. As we did in the case of the approach 
based on the corrector equation, a first step towards the analysis of this numerical 
method is to quantify the corresponding qualitative result, namely here Kipnis- 
Varadhan's convergence theorem. Compared to the deterministic approach based 
on the approximate corrector equation, the advantage of the present approach is 
that its convergence rate and computational costs are dimension-independent. As 
we shall also see, as opposed to the approach based on the corrector equation, the 
environment only needs to be generated along the trajectory of the random walker, 
so that much less information has to be stored during the calculation. This may be 
quite an important feature of the Monte Carlo method in view of the discussion of 
[GTPOl Section 4.3]. 

We consider the discrete elliptic operator —V* ■ AV, where V*- and V are the 
discrete backward divergence and forward gradient, respectively. For all x G Z d , 
A(x) is the diagonal matrix whose entries are the conductances w x . x + ei of the edges 
(x,x + ej) starting at x, where (e»)ie{l,...,d} denotes the canonical basis of Z d . Let 
B denote the set of edges of Z d . We call the family of conductances uj = (w e ) e eB 
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the environment. The environment w is random, and we write P for its distribution 
(with corresponding expectation E). We make the following assumptions : 

(HI) the measure P is invariant under translations, 
(H2) the conductances are i. i. dQ, 

(H3) there exists < a < (3 such that a ^ ui e ^ (3 almost surely. 

Under these conditions, standard homogenization results ensure that there exists 
some deterministic symmetric matrix Ah om such that the solution operator of the 
deterministic continuous differential operator —V • AhomV describes the large scale 
behavior of the solution operator of the random discrete differential operator —V* • 
AV almost surely (for this statement, (H2) can in fact be replaced by the weaker 
assumption that the measure P is ergodic with respect to the group of translations, 
see [Ku83] L 

The operator —V* • AV is the infinitesimal generator of a stochastic process 
(X(t))teR + which can be defined as follows. Given an environment w, it is the 
Markov process whose jump rate from a site x £ Z d to a neighbouring site y is 
given by uj x ,y We write P" for the law of this process starting from x £ Z d . 

It is proved in [KV86] that under the averaged measure PPg , the rescaled pro- 
cess y/eX(e~ 1 t) converges in law, as e tends to 0, to a Brownian motion whose 
infinitesimal generator is —V • AhomV, or in other words, a Brownian motion with 
covariance matrix 2Ah om (see also [AKS82, Kii83j IKo85] for prior results). We 
will use this fact to construct computable approximations of A\ lom . As proved in 
|DFGW8"9] . this invariance principle holds as soon as (HI) is true, (H2) is replaced 
by the ergodicity of the measure P, and (H3) by the integrability of the conduc- 
tances. Under the assumptions (H1-H3), SS04] strengthens this result in another 
direction, showing that for almost every environment, s/eX(£ t) converges in law 
under Pq to a Brownian motion with covariance matrix 2Ay lom . This has been it- 
self extended to environments which do not satisfy the uniform ellipticity condition 
(H3), see [BB071 IMP071 IBP071 lMa08l IBDIO] . 

Let (Y(t)) t& N denote the sequence of consecutive sites visited by the random walk 
(X(t))teR + (note that the "times" are different in nature for X(t) and Y(t)). This 
sequence is itself a Markov chain that satisfies for any two neighbours x,y £ Z d : 



P£[F(l)=y] = 



Pu>{x) 

where Puj(x) = J2\z\=i ^x^+z- We simply write p{u>) for p w (0). Let us introduce a 
"tilted" version of the law P on the environments, that we write P and define by 

(1.1) dPH = ?M dPH . 

The reason why this measure is natural to consider is that it makes the environment 
seen from the position of the random walk Y a stationary process (see (|3.2p for a 
definition of this process). 

Interpolating between two integers by a straight line, we can think of Y as a 
continuous function on M + . With this in mind, it is also true that there exists a 
matrix A^ 1 ^ such that, as e tends to 0, the rescaled process y / e : ^(£ _1 converges 
in law under PPq to a Brownian motion with covariance matrix 2^^. Moreover, 
A hom and ^hom are related by (see [DFGW891 Theorem 4.5 (ii)]) : 

(1.2) A hom = E\p]A^ 



\ disc 
horn " 



(H2) obviously implies (HI) in the present form. Yet for most qualitative (and some quanti- 
tative) results (H2) can be weakened and may not imply (HI) any longer. 
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Given that the numerical simulation of Y saves some operations compared to the 
simulation of X (there is no waiting time to compute, and the running time is equal 
to the number of steps), we will focus on approximating Ajjjf^. More precisely, we 
fix once and for all some £ € K d with |£| = 1, and define 

(1.3) o? = r 1 EESf[(e.y(t)) 2 ], a 2 = 2^.Af^. 

It follows from results of |KV86j (or [DFGW891 Theorem 2.1]) that of tends to a 2 
as t tends to infinity. Our first contribution is to give a quantitative estimate of 
this convergence. In particular we shall show that, with i. i. d. coefficients and up 
to a logarithmic correction in dimension 2, the difference between a 2 and a 2 is of 
order 1/t. 

We now describe a Monte-Carlo method to approximate of. Using the definition 
of the tilted measure (jl.ip , one can see that 

, _ EE£[(£ • Y(t)) 2 ] EE£[p( W )(£ • Y{t)) 2 ] 



(1.4) 



t m\p] 



Assuming that we have easier access to the measure P than to the tilted P, we 
prefer to base our Monte-Carlo procedure on the r. h. s. of the second identity in 
(|1.4p . Let yw, Y^ 2 \ ... be independent random walks evolving in the environments 
ujM 7 a/ 2 ) j . . . respectively. We write Pq for their joint distribution, all random walks 
starting from 0, where uj stands for (u^' , lo^ 2 \ . . .) . The family of environments ZJ is 
itself random, and we let P® be the product distribution with marginal P. In other 
words, under P®, the environments wW,a/ 2 ), . . . are independent and distributed 
according to P. Our computable approximation of a 2 is defined by 

. p(q,W)(g ■ YW(t)) 2 + ■ ■ ■ + p(c»)(g ■ Y^(t)) 2 

1 J t(p(w(i)) + ...+p(wW)) 

The following step in the analysis is to quantify the random fluctuations of A n {t) 
in terms of n — the number of random walks considered in the empirical average 
to approximate erf — and t. We shall prove a large deviation result which ensures 
that the P®Pq -probability that the difference between A n (t) and of exceeds 1/t is 
exponentially small in the ratio n/t 2 . 

The rest of this article is organized as follows. In Section [5J which can be read 
independently of the rest of the paper, we consider a general discrete or continuous- 
time reversible Markov process. Kipnis-Varadhan's theorem (and its subsequent 
development due to [DFGW89 ) gives conditions for additive functionals of this 
process to satisfy an invariance principle. We show that, under additional conditions 
written in terms of a spectral measure, the statement can be made quantitative. 
More precisely, Kipnis-Varadhan's theorem relies on writing the additive functional 
under consideration as the sum of a martingale plus a remainder. This remainder, 
after suitable normalization, is shown to converge to in L 2 . Under our additional 
assumptions, we give explicit bounds on the rate of decay. In Section^ we make use 
of this result, in the context of the approximation of homogenized coefficients, to 
estimate the systematic error \a 2 — a 2 \ . The central achievement of this section is to 
prove that the relevant spectral measure satisfies the conditions of our quantitative 
version of Kipnis-Varadhan's theorem. Section fj] is dedicated to the estimate of 
the random fluctuations. These are controlled through large deviations estimates. 
Relying on these results, we give in Section [5] a complete error analysis of the Monte- 
Carlo method to approximate the homogenized matrix Afj^, which we illustrate 
by numerical tests. 

Let us quickly discuss the sharpness of these results. If A was a periodic matrix 
(or even a constant matrix) the systematic error would also be of order 1 /t (without 
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logarithmic correction for d — 2), and the fluctutations would decay exponentially 
fast in the ratio n/t 2 as well. This shows that our analysis is optimal (the additional 
logarithm seems unavoidable for d — 2, as discussed in the introduction of |GO10aj h 

Let us also point out that although the results of this paper are proved under 
assumptions (H1)-(H3), the assumption (H2) on the statistics of uj is only used to 
obtain the variance estimate of [G OlOal Lemma 2.3]. In particular, (H2) can be 
weakened as follows: 

• the distribution of W( Z)Z -|_ ej ) may m addition depend on e^, 

• independence can be replaced by finite correlation length Cl > 0, that is 
for all e, e' S B, w e and u) e > are independent if je — e'j ^ Cl, 

• independence can be replaced by mixing in the sense of Dobrushin and 
Shlosman — we refer the reader to work in progress by Otto and the first 
author for this issue |G011j . 

Notation. So far we have already introduced the probability measures Pg (dis- 
tribution of Y), Pq 7 (distribution of Y {2 \ . . .), P (i.i.d. distribution for lu = 
(w e ) ee s), P (tilted measure defined in (jl.llO and P® (product distribution of uJ with 
marginal P). It will be convenient to define P® the product distribution of 57 with 
marginal P. For convenience, we write Po as a short-hand notation for PPq , Pq for 
PPft, Pf for P®Pff, and Pjf for P^Py 7 . The corresponding expectations are written 
accordingly, replacing "P" by "E" with the appropriate typography. Finally, we 
write | • | for the Euclidian norm of M. d . 

2. Quantitative version of Kipnis-Varadhan's theorem 

Kipnis-Varadhan's theorem |KV86j concerns additive functionals of reversible 
Markov processes. It gives conditions for such additive functionals to satisfy an 
invariance principle. The proof of the result relies on a decomposition of the additive 
functional as the sum of a martingale term plus a remainder term, the latter being 
shown to be negligible. In this section, which can be read independently of the rest 
of the paper, we give conditions that enable to obtain some quantitative bounds on 
this remainder term. 

We consider discrete and continuous times simultaneously. Let (^t)t^o be a 
Markov process defined on some measurable state space H (here, t ^ stands either 
for t € N or for t 6 R+). We denote by P x the distribution of the process started 
from x <G H, and by E x the associated expectation. We assume that this Markov 
process is reversible and ergodic with respect to some probability measure v. We 
write P v for the law of the process started from the distribution and E v for the 
associated expectation. 

To the Markov process is naturally associated a semi-group (Pt)t>o defined, for 
any / G L 2 (f), by 

P t f(x) = E x [f(Vt)]- 

Each P t is a self-adjoint contraction of ~h 2 (v). In the continuous-time case, we 
assume further that the semi-group is strongly continuous, that is to say, that P t f 
converges to / in L 2 (^) as t tends to 0, for any / S ~h 2 {v). We let L be the h 2 (v)- 
infinitesimal generator of the semi-group. It is self-adjoint in L 2 (z/), and we fix the 
sign convention so that it is a positive operator (i.e., P t — e~ tL ). 

Note that in general, one can see using spectral analysis that there exists a 
projection P such that P t f converges to Pf as t tends to 0, t > 0. Changing L 2 (i/) 
to the image of the projection P, and Po for P, one recovers a strongly continuous 
semigroup of contractions, and one can still carry the analysis below replacing L 2 (^) 
by the image of P when necessary. 
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In discrete time, we set L = Id — Pi. Again, L is a positive self-adjoint operator 
on L 2 (V). Note that we slightly depart from the custom of defining the generator 
as Pi in order to match more closely the continuous time situation. 

We denote by (•,•) the scalar product in L 2 (z/). For any function / e L 2 (^) 
we define the spectral measure of L projected on the function / as the measure e/ 
on R + that satisfies, for any bounded continuous \& : R+ — > R, the relation 

(2.1) </,*(£)/} = J *(A)de/(A). 
The Dirichlet form associated to L is given by 

(2.2) \\f\\j= J \de f (\). 

We denote by H 1 the completion of the space {/ S L 2 (f) : ||/||i < +00} with 
respect to this || • ||i norm, taken modulo functions of zero || • ||i norm. This turns 
(H , || • ||i) into a Hilbert space, and we let H" 1 denote its dual. One can identify 
H^ 1 with the completion of the space {/ £ L 2 (^) : ||/||-i < +00} with respect to 
the norm |j • ||_i defined by 

ll/ll-i= /A" 1 d e/ (A). 
Indeed, for all / £ L 2 (^), the linear form 

/ (L»nff\|| - Hi) R 
I <t> h> </,<*> 

has norm ||/||_i, and thus defines an element of H^ 1 (with norm ||/||_i) iff ||/||-i 
is finite. The notion of spectral measure introduced in (|2.1j) for functions of H. 2 (i/) 
can be extended to elements of Indeed, let ^ : R+ — > R be a continuous 

function such that ^(A) = 0(A _1 ) as A — > +00. One can check that the map 

/ (iL»n#-\|HI-i) tf 1 
I / ^ *(£)/ 

extends to a bounded linear map on H . One can then define the spectral measure 
of L projected on the function / as the measure e/ such that for any continuous ^ 
with *(A) = 0(A _1 ), (|2~Tj) holds. With a slight abuse of notation, for all / e H^ 1 
and g <E H , we write (/,<?) for the — P 1 duality product between / and g. 
For any / 6 if -1 , we define (Zf(t)) t ^o as 

/■* t-l 
(2-3) Z/(t) =//(%) ds or Z/(t) = 

according to whether we consider the continuous or the discrete time cases. In 
the continuous case, the meaning of (12.31) is unclear a priori. Yet it is proved in 
[DFGW891 Lemma 2.4] that for any t > the map 

f ^Mnr 1 -> l 2 (p„) 
I / -> Z/(t) 

can be extended by continuity to a bounded linear map on H , and moreover, 
that (|2.3[) coincides with the usual integral as soon as / G L 1 ^). The following 
theorem is due to [DFGW89 , building on previous work of |KV86) . 
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Theorem 2.1. (i) For all f € H^ 1 , there exists (Mt)t^o, (ft)t^o such that Zf(t) 
defined in (|2.3p satisfies the identity Zf(t) = M t + ft, where (M t ) is a square- 
integrable martingale with stationary increments under P v ( and the natural filtra- 
tion), and (ft) is such that : 

(2.4) t - 1 E4(tt)*]——+0. 

As a consequence, t~ x / 2 Z f(i) converges in law under P v to a Gaussian random 
variable of variance o~ 2 {f) < +oo as t goes to infinity, and 

(2.5) t -i Eu[{Zf (t)?]—-^o 2 {f). 

t— y+oo 

(ii) If, moreover, f £ L 1 (^) and, for some t > 0, sup ^ t ^( |^/(i)l * s ^ n L 2 (f), then 
the process t i— > \feZ f{e t) converges in law under P v to a Brownian motion of 
variance o 2 (f) as e goes to 0. 

Remarks. The additional conditions appearing in statement (ii) are automatically 
satisfied in discrete time, due to the fact that C h 2 (v) in this case. In the 

continuous-time setting and when / e L 1 (^), the process 1 1-> Zf(t) is almost surely 
continuous, and sup ^ ts jj |^/(*)| is indeed a well-defined random variable. 

Under some additional information on the spectral measure of /, we can estimate 
the rates of convergence in the limits (|2.4I) and (I2.5[) . For any 7 > 1 and q ^ 0, we 
say that the spectral exponents of a function / e H^ 1 are at least (7, —q) if 

(2.6) [" 'de f (\)=0(^]n«( f x- 1 )) (/x^O). 

JQ 

Note that the phrasing is consistent, since if (7', —q') < (7, —q) for the lexicograph- 
ical order, and if the spectral exponents of / are at least (7, — q), then they are 
at least (7',— 9')- I n |Mollj . it was found more convenient to consider, instead of 
(I2.6[) . a condition of the following form : 

(2.7) ^ \- x de/(A) = 0{ff- x \a 9 {jx- x )) (/i^O), 
Jo 

One can easily check that conditions (|2.6[) and (|2.7I) are equivalent. Indeed, on the 
one hand, one has the obvious inequality 

de/(A) <p A" 1 de/(A), 
Jo 

which shows that (|2.7[) implies (|2.6[) . On the other hand, one may perform a kind 
of integration by parts, use Fubini's theorem : 

A" 1 de/(A) = / / S- 2 dS de/(A) 

<5A/i 

2 ' 



S- 2 / de/(A) dJ, 
5=0 J\=o 

and obtain the converse implication by examining separately the integration over 
S in [0,/i) and in [fi, +00). 



For all 7 > 1 and g ^ 0, we set 
(2.8) V 7 , 9 W - 

The quantitative version of Theorem 12.11 is as follows 



i" 1 ^ 1 ^) if 7 = 2, 
i" 1 if 7 >2. 
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Theorem 2.2. If the spectral exponents of f £ H~ x are at least (7, —q), then the 
decomposition Zf(t) = M t + f % of Theorem \2.1\ holds with the additional property 
that 

t-'EMt) 2 } =0(V 7 , g (t)) (t->+oo). 

Moreover, 

Proof. In the continuous-time setting, the argument for the first estimate is very 
similar to the one of [Molll Proposition 8.2], and we do not repeat the details here. 
It is based on the observation that 

(2-9) i^[(6) 2 ] = 2 m^ de/ (A). 

One needs to take into account the possible logarithmic terms that appear in (I2.7[) 
and which are not considered in [Mollj . Some care is also needed because we do 
not assume that / G L 2 (V). Yet one can easily replace the bound involving the 
L 2 (z/) norm of / by its H^ 1 norm. The second part of the statement is given by 
[Molll Proposition 8.3]. 

We now turn to the discrete time setting. In this context, identity (|2.9p should 
be replaced by 

1, r„*»i_n H-a-A)' 



^MtV] = 2 J de/(A). 

By definition, L = Id — Pi, where Pi is the semi-group at time 1. Hence the 
spectrum of L is contained in [0, 2]. One can then follow the same computations as 
before to prove the first part of Theorem 12.21 

Somewhat surprisingly, the second part of the statement requires additional at- 
tention in the discrete time setting. Indeed, in the continuous case, the argument 
of [Molll Proposition 8.3] (which already appears in [DFGW89 ) is that Zj{t) and 
are orthogonal in L 2 (P„), a fact obtained using the invariance under time sym- 
metry. This orthogonality is only approximately valid in the discrete-time setting. 
Indeed, let us recall that Zf(t) is given by (|2.3[) . while £t is obtained as the limit in 
L 2 (P„) of 

-u e (r) t ) + u £ (r] ), 

where u s = (e + L)^ 1 f. Using time symmetry, what we obtain is that £t is orthog- 
onal to (Zf(t) + f{r) t ))- As a consequence, the cross-product E v [Zf(t)^ t ], which is 
equal to in the proof of [Molll Proposition 8.3], is in the present case equal to 
— E„[f(r)t)£t\- Yet spectral analysis ensures that this term is equal to 

l-(l-A) de/(A) = Q(1) (t^+00), 
A 

which is what we need to obtain the second claim of the theorem. □ 



3. The systematic error 

We now come back to the analysis of the Monte-Carlo approximation of the 
homogenized coefficients within assumptions (H1)-(H3). The aim of this section is 
to estimate the difference between of and the quantity a 2 we wish to approximate 
(both being defined in (|1.3I) V This difference, that wc refer to as the systematic 
error after [GOlOaj . is shown to be of order 1/t as t tends to infinity, up to a 
logarithmic correction in dimension 2. 
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(3.1) 



Theorem 3.1. Under assumptions (Hl)-(HS), there exists q such that, as t 
tends to infinity, 

0(t- l W{t)) ifd = 2, 

olt- 1 ) ifd>2. 

Theorem 13. II is a discrete-time version of |Molll Corollary 2.6]. Its proof makes 
use of an auxiliary process that we now introduce. 

Let (0 x ) X £Z d be the translation group that acts on the set of environments as fol- 
lows: for any pair of neigbhours y,z € Z rf , (6 X cj) v , z = 0J x +y,x+z- The environment 
viewed by the particle is the process defined by 

(3.2) w(t) = 6 Y(t) u. 

One can check that (ui(t))t£N is a Markov chain, whose generator is given by 

1 x 



(3.3) 



wo,z{f(6, u)-f{u)), 



so that Eft[/(w(l))] = (I - C)f{uj). Moreover, the measure P defined in (fTTTj) is 
reversible and ergodic for this process IDFGW89] Lemma 4.3 (i)] . As a consequence, 
the operator C is (positive and) self-adjoint in L 2 (P). 

The proof of Theorem 13 . 1 1 relies on spectral analysis. For any function / E L 2 (¥), 
let e/ be the spectral measure of C projected on the function /. This measure is 
such that, for any positive continuous function ^ : [0, +oo) — > R+, one has 

E[/ *(£)/] - J *(A) d e/ (A). 

For any 7 > 1 and q ^ 0, we recall that we say that the spectral exponents of a 
function / are at least (7, —q) if (|2.6p holds. 
Let us define the local drift in direction £ as 

a(w) = E5fK-y(i)] = -L 



(3.4) 



As we shall prove at the end of this section, we have the following bounds on the 
spectral exponents of 0. 

Proposition 3.2. Under assumptions (Hl)-(HS), there exists q ^ such that the 
spectral exponents of the function d are at least 

(2,-q) ifd=2, 

(d/2 + 1,0) i/3<d<5, 

(4,-1) ifd=6, 

(4,0) ifd^l. 

Let us see how this result implies Theorem 13. II In order to do so, we also need 
the following information, that is a consequence of Proposition [XU 

Corollary 3.3. Let 

(3.5) t> t (u) = Egf[9(«(t))] 

be the image of d by the semi-group at time t associated with the Markov chain 
(uj(t))t£®- There exists q ^ such that 

0(t- 2 ln q (t)) ifd = 2, 
) t - W 2+i)\ 

0(t- 4 ln(i)j ifd = 6, 
oh' 4 ) ifd>7. 



E[(0 t 
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Proof. This result is the discrete-time analog of |GM1Q[ Corollary 1] . It is obtained 
the same way, noting that 

%>*) 2 ]= |(l-A) 2t de f (\), 

and that the support of the measure e/ is contained in [0, 2]. □ 
We are now in position to prove Theorem 13. II 



Proof of Theorem \3.1l The proof has the same structure as for the continuous-time 
case of [Molll Proposition 8.4]. Note that |DFGW891 Theorem 2.1] ensures that 

(3.6) lim a 2 ( = f) lim ■ Y(t)) 2 ] = a 2 . 

The starting point is the observation that, under Pq, the process defined by 

t—i 

(3.7) N t = ^Y(t)-^D(oj(s)) 

s=0 

is a square-integrable martingale with stationary increments. On the one hand, 
following (I2.3[) . we denote by Z$(t) the sum appearing in the r. h. s. of (13.71) . From 
Proposition 13.21 and Theorem 12.21 we learn that there exist a and q ^ such that 



(3.8) ta 2 -E [(Z d (t)) 2 



0(ln q (t)) itd =2, 
0(1) if d > 2. 



On the other hand, since N t is a martingale with stationary increments, 
(3.9) M(N t ) 2 } = tE Q [(N 1 ) 2 }. 

As in the proof of Theorem 12.21 in the discrete time case, we then use that £ • Y(t) 

is orthogonal to (Z$(t) + Z)(w(f))) to turn (|3.7[) into 

(3.10) 

*- x Eo[(iV t ) a ] = r x E [(£ ■ Y(t)) 2 } + t^EoKZ^t)) 2 } + 2^ 1 E [5(c(t))(e • Y {£))]. 

We already control the 1. h. s. and the second term of the r. h. s. of (13.10[) . In 
order to quantify the convergence of £ -1 Eo[(£ • Y(t)) 2 ] it remains to control the last 
term. In particular, provided we show that 



(3.11) Eo[3(w(t))(£-r(t))] = 



0(ln q (t)) if d = 2, 
0(1) if d> 2, 



([3T0]) . (HH) , dSH) , and (JSH) imply first that a 2 = E [(7Vi) 2 ] - a 2 , and then the 
desired quantitative estimate (|3.ip . We now turn to (|3.11l) and write 

t-i 

E mu;(tm-Y(t))} = ^Eo[D( W (t))(€.(F( a + l)-y( fl ))] 

t-1 

= 5^Eo[9 t _ a _i(w(s + 1))(C • (^(fl + 1) - 

where we have used the Markov property at time s + 1, together with the definition 
(I3.5[) of t _ s _i. Using Cauchy-Schwarz inequality and the stationarity of the process 
(w(t))t e N under Eo, this sum is bounded by 

|C| 2 EE[(Wi) 2 ] 1/2 . 

Esimate p. lip then follows from Corollary 13.31 This concludes the proof of the 
theorem. □ 
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We finally turn to the proof of Proposition 13.21 which is a discrete-time coun- 
terpart of GMlQj Theorem 5]. In jQMlQl Theorem 5] however, we had proved in 
addition that the spectral exponents are at least (d/2 — 2, 0), which is sharper than 
the exponents of Proposition 13.21 for d > 10. In particular for d > 10 the bounds 
of [GM101 Theorem 5] follow from results of Moll , whose adaptation to the dis- 
crete time setting is not straightforward. As shown above, the present statement 
is sufficient to prove the optimal scaling of the systematic error, and we do not 
investigate further this issue. The proof of Proposition 13.21 is rather involved and 
one may wonder whether this is worth the effort in terms of the application we have 
in mind — namely Theorem 13.11 In order to obtain the optimal convergence rate 
in Theorem 13.11 we need the spectral exponents to be larger than (2,0). Proving 
that the exponents are at least (2,0) is rather direct using results of |GO10aj (see 
the first three steps of the proof of Proposition 13. 2[) . Yet proving that they are 
larger than (2,0) for d > 2 is as involved as proving Proposition 13.21 itself. This is 
the reason why we display the complete proof of Proposition 13.21 — although the 
precise values of the spectral exponents are not that important in the context of 
this paper. 

Proof of Proposition VS. "A The present proof is a direct version of the proof of GM10 , 
Theorem 5]. In particular, the proof of [GM101 Theorem 5] relies on a nontrivial 
(weaker) estimate of the spectral exponents obtained in [GOlObl using a covari- 
ance estimate. Yet if one wants to extend these results to more general statistics 
of the conductivity function — for instance to mixing coefficients in the sense of 
Dobrushin and Shlosman — one has to give up the covariance estimate (which we 
are not able to prove any longer, see |GQllj ). With this in mind we only rely 
on the variance estimate of [GOlOa, Lemma 2.3]. Our strategy is similar to the 
strategy used for (continuous) elliptic equations in |GQllj to prove the estimate 
of the systematic error. In particular we directly focus on the spectral exponents 
rather than on some other quantity like the systematic error itself. Yet we will go 
slightly further. Starting point is the inequality: 

(3-12) J de,(X)<T- 4 J^ +A)4 de (A), 

which follows from the fact that for A ^ T _1 , (t-i + x)a > 1 (here and below 
< and > stand respectively for ^ and up to multiplicative constants). The 
variable T -1 for T large plays the role of fi in (|2.6p . In what follows we make the 
standard identification between stationary functions (z,uj) i-> f(z,(j) of both the 
space variable z e 7L d and the environment u and their translated versions at 
u> i y /(0, z u) depending on the environment only. We define <f>T as the unique 
stationary solution to 

(3.13) T-Vt(x) - ~4t V * • A(x)W4>t = 4t V * ' A (*)Z, 

whose existence and uniqueness follow from Lax-Milgram's theorem in L 2 (P) using 
the identification between the stationary function <pT and its version defined on 
the environment only (see a similar argument of Kii83]). In particular, with the 
notation d = V* • A(x)£, 

4>t = (T" 1 + 

where C is the operator defined in (|3.3p . and the spectral theorem ensures that 

f°° 1 

E(4) = EMT- 1 + C)- 2 *) = y o (T _ 1 + A)2 de a (A) 
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where e is the spectral measure of L projected on the drift 0. We also let ipr be 
the unique stationary solution to 



(3.14) 



T"Vt(z) 



1 



V* • A(x)Vip T (x) = <fo{x), 



whose existence and uniqueness also follows from Lax-Milgram's theorem in the 
probability space as well. This time, 

i'r = {t~ x + c)-\ 

and the spectral theorem yields 

E(V4) = E^T- 1 + £)- 4 0) = 



1 



(T-i + A) 



jde (A) 



From now on, we shall use the shorthand notation (u) := E(u) and var [u] 
((u — (w)) 2 ) for all u G L 2 (P). In particular the identity above turns into 

1 



(3.15) 



deo(A) = var [ipq 



(T- 1 +A) 4 

since (ipr) = fijj / i>TV dP = / 4>tP dP = using equations ({3.140 and (|3.13p . 

The rest of the proof, which is dedicated to the estimate of var [ipr] , is divided 
in six steps. Starting point is the application of the variance estimate of [GOlOal 
Lemma 2.3] to ipT, which requires to estimate the susceptibility of ipT with respect 
to the random coefficients. In view of (|3.14j) it is not surprising that we will have 
to estimate not only the susceptibility of ipx but also of 4>t and of some Green's 
function with respect to the random coefficients. In the first step, we establish 
the susceptibility estimate for the Green's function. In Step 2 we turn to the 
susceptibility estimate for the approximate corrector <^j\ We then show in Step 3 
that, relying on {GOlOaj , this implies that the spectral exponents are at least 



(3.16) 



(2, -q) 
(2,0) 



for 
for 



d = 2, 
d>2. 



Following step is to estimate the susceptibility of -0T- In Step 5 we show that this, 
combined with the suboptimal estimates 
(3.17) 



<IW>t| 2 ) 



< 



Tin 9 T 
T 



for 
for 



d = 2, 
d>2, 



and (ipx) 



< 



T 2 In" T 
T 2 



for 
for 



d=2, 
d>2 



obtained using Q3. 160 . allow us to improve the spectral exponents to 
(3.18) 



(2,-9) 


for 


d = 2, 


(5/2,0) 


for 


d = 3, 


(3,-1) 


for 


d = 4, 


(3,0) 


for 


d > 4. 



In the last step, we quickly argue that in turn these spectral exponents yield the 
optimal and suboptimal estimates 
(3.19) 



(IVVtI 2 ) 



< 



Tin 9 T 

VT 

InT 
1 



for 
for 
for 
for 



d = 2, 
d = 3, 
d = 4, 
d > 4, 



and ("0t) 



< 



T 2 ln 9 T for (2 = 2, 

T 3 / 2 for (2=3, 

TlnT for (2=4, 

T for d > 4, 
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which finally bootstrap (|3.18[) to the desired estimates of the spectral exponents, 
and consequently yield the following optimal estimate of var [?/>t] : 



var [ipr] 



< 



T 2 ln"T 


for 


d = 2, 


T 3/2 


for 


d = 3, 


T 


for 


d = 4, 


VT 


for 


d = 5, 


InT 


for 


d = 6, 


1 


for 


d>6. 



We do not need Step 6 to prove Theorem 13.11 Yet this step is interesting in itself 
since it may be the starting point of a subtle "multi-level" induction to obtain the 
sharp spectral exponents in any dimension. 

Step 1 . Susceptility of the Green's function. 

For all y 6 Z d we define the Green's function Gt(-, y) with singularity at y as the 
unique solution in L 2 (Z d ) to the equation 



(3.20) 



T^p^Grix, y) - VJ • A(x)V x G T (x, y) = 5(x - y), 



using Lax-Milgram's theorem. We shall prove for all e = (z, z') € B, z € Z d , 

z' = z + e», 

(3.21) ^(x,y) = -T- 1 ( K G T (z,y)GT(x,z) + G T (z',y)GT(x,z')) 

-V Zi G T (x,z)X7 Zi G T (z,y), 



and 
(3.22) 

where gx : Z d 
(3.23) 

for d > 2, and 
(3.24) 



sup\V Zi G T (z,y)\ < \V Zi G T {z,y)\+T- l g T (y-zl 
sup\V Zt G T (y,z)\ < \V Zi G T (y,z)\+T- l g T (y-z), 



satisfies for some constant c > (depending on a, (3, d) 
g T {x) = (l + | a; |) 2 - d cxp(- c ^) 



9t(x) 



ln (iTld ) 



/ F \ 

exp ^ ~ c Vt* 



for d = 2. We define the elliptic operator Lt as 

(L T u)(a;) = uj^^T^uix) + ^(s,*') ("(a;) - u(x')), 

x' ,\x — x'\ — 1 x',|a;— a;'|=l 

so that f|3.20p takes the form 

(3.25) (L T G T (;y))(x) = 5{x-y). 

Formally differentiating this equation with respect to u e yields 
'dGi 



L 



t(^(-,2/)) (x) +T- 1 G T (x,y)(S(x - z) + S(x - z')) 



+ (G T (z, y) - G T (z', y))S(x - z) + (G T (z', y) - G T (z', y))S(x - z') = 0. 
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Using (|3.25[) this identity turns into 

Lt (^(-, y) + T-^Griz, y)G T (; z) + G T (z' \ y)G T (; *')) 

+ V Zi G T {;z)V Zi G T {z,y)\(x) = 0. 

Formally erasing the operator Lt on the 1. h. s. then yields the desired identity 
(|3.2ip . To turn this into a rigorous argument, we may proceed as in |GO10a[ Proof 
of Lemma 2.5], first consider finite differences instead of derivative w. r. t. ui e , use 
that Lt is bijective on L 2 (Z d ), and then pass to the limit. We leave the details to 
the reader and directly turn to (|3.22p . 
From (|3.21[) we infer that 



dV Zi G T (z,y) 



-V Zi V Zi G T (z,z)V Zi G T (z,y) 

- T- l {G T (z, y)V Zt G T (z, z) + G T (z', y)V Zi G T {z, z')) . 



Using then the uniform pointwise estimate |VGt| 1 (see [GOlOal Corollary 2.3], 
whose proof holds as well in the present case with a non-constant zero order term 
since 2 d aT~ 1 ^ T~ 1 p u (x) ^ 2 d f3T~ 1 ), the uniform pointwise estimate of the 
Green's function from |G110, Lemma 3], and considering this identity as an ODE 
for V Zi GT(z,y) in function of u e , we obtain (I3.22p . 

Step 2. Susceptibility of the approximate corrector (f>T- 

In this step we shall prove that for e = (z, z') € B, z £ Z d and z' = z + e,, 



(3.26) 



duj. 



■(x) = -(Vi4> T (z)+ti)V Zi G T (x,z) 

-T- 1 Mz)(Gt(x,z)+G t (x,z')), 



(3.27) sap\4>r(x)\ £ \M*)\ 

+ (|Vi^r(2)| + l)(|V^G T (x, z)\ + T-^ 2 g T (x - z)) , 



(3.28) sup 



-Or) 



8uj l 

and for all n e N, 

d{(t> T {x) n+1 



< (|V<£r(*)l + l)(\V Zi G T (x,z)\ + T-V 2 g T (x - z)), 



(3.29) sup 



8uJ e 



< (|V<M*)I + l)(\V Zi G T (x, z)\ + T-V 2 g T {x - z)) 

x (\4>t(x)\ + (|Vi<M*)| + l)(\V Zi G T (z,x)\ + T-^ 2 g T (x - z))". 
As for the Green's function, we rewrite the defining equation for 4>t as 
(3.30) (L T G T (;y)){x) - V* ■ A(x)£ = 0. 

Formally differentiating (|3.30[) w. r. t. u> e yields 

Lt ^7 {x) ~ ^^ T ^ + - z ) - 5 i x - 

+ rVr(z)(<S(i -z) + 8(x - z')) = 0, 
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which using p.25[) turns into 



duj e 



<yiM*) + Zi)(G T (;z)-G T (;z')) 



+ T- 1 c/ )T (z)(G T (-,z) + G T (-,z , ))yx) = 0. 



This (formally) shows p.26[) . To turn this into a rigorous argument, we may use 
the Green representation formula 



cj> T {x) = I G T (x,y)V* ■ A(y)£ dy, 



where L d dy denotes the sum over all y €E Z d , and proceed as in |GO10a[ Proof of 
Lemma 2.4]. 

We now turn to ([3728]) . This estimate follows from (pT26]> . p~22|) . and the fol- 
lowing two facts: 



(3.31) 
and 
(3.32) 



\<h\ < Vt 



sup|Vi^r(z)| < \Vi(j> T {z)\ + 1. 
Starting point to prove (|3.3ip is the Green representation formula in the form of 



I Ms) | 



G T (x,y)V* • A(y)i dy 



< 



f V y G T (x,y) ■ A(y)£ dy 
V y G T (0,y)| dy, 



from which we deduce the claim by a dyadic decomposition of space combined 
with Cauchy-Schwarz' inequality and jGOlOai Lemma 2.9] (a similar calculation is 
detailed in jGUOl Proof of Lemma 4]). For (I3.32j) . we first note that p. 261) implies 



= -(ViMz) + £i)(V, i G T (* / , z) - V 2i G T (z, z)) 

+ T~ x ch{z) (V„G r (*, z) + V z< G T (z, z')) , 



which — seen as an ODE w. r. t. cj e — yields the claim using the uniform bound 
|VG T | < 1 of IGOlOal Corollary 2.3] and (pOT) . 

Estimate (|3.27p is a direct consequence of p.28p . whereas p. 291) follows from 
the Leibniz' rule combined with pUSJ), (j3~?f| . and p~2"g]) . 

Step 3. Proof of (|3~16)) . 

The estimates p,16[) of the spectral exponents follow from the more general esti- 
mates: for all q > there exists 7(g) > such that 

ln 7(<?) T for d= 2, 
1 for d > 2, 



(3.33) 



<I<H ? ) 



< 



combined with the fact that 



pi pC 

/ de (A) <T- 2 f 
Jo Jo 



1 

(T-i + xy 



de D (A) = T- 2 <4) 
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The proof of (j3.33[) is an adaptation of GOlOa, Proof of Proposition 2.1] which 
already covers the case of a constant coefficient in the zero order term of Lt, that 
is for T~ x (\)t instead of T~ 1 p UJ 4>T (no randomness in the zero order term). 

The first step to apply the variance estimate of [GOlOal Lemma 2.3] is to show 
that cf>T is measurable with respect to cylindrical topology associated with the 
random variables. This can be proved exactly as in [G OlOal Lemma 2.6]. The 
estimates (|3~2"2"|) . (|3~27)l . (pOSl) . and (|3T2"9")) of Steps 1 and 2 are like in the auxiliary 
lemmas of [GOlOa provided we replace the terms \V Zi Gx(x, z)\ in [GOlOal Lem- 
mas 2.4 & 2.5] by \V Zi G T (x, z)\ + T^ 1 l 2 g T {x - z). A close look at the proof of 
[GOlOal Proposition 2.1] shows that these terms \V Zi GT(x, z)\ are either estimated 
by the Green's function itself (in which case the additional term T~ x / 2 gT(x — z) is 
of higher order), or they are controlled on dyadic annuli by the Meyers' estimate 
[GOlOal Lemma 2.9]. This lemma shows in particular that there exists p > 2 such 
that for all p ^ q ^ 2, k > and R » 1, 

\V z G T (x,z)\ q dz < R d (R 1 ~ d ) q mm{l,VTR- 1 } k . 

I R^\x-z\<2R 

By the properties (13.23P for d > 2 and (j3.24[) for d = 2 of the function qt , it is easy 
to see that for d > 2 

{T- l l 2 g T (x- z)) q dz < R d {R 1 - d ) q mm{\,VTR- 1 } k 

I R^\x-z\<2R 

as well, whereas for d = 2 

(T-^ 2 g T {x~ z)) q dz < R 2 (R- 1 ) q ln q Tmm{l,VfR- 1 } k . 

I R^\x-z\<2R 

Hence the proof of GOlOa, Proposition 2.1] adapts mutadis mutandis to the present 
case (with possibly an additional logarithmic correction for d = 2), and we have 
(15331 . 

Step 4. Susceptibility of ipT- 

In this step we shall prove that for e = (z, z'), z <E % d and z' = z + e^, 
(3.34) ^Z(x) = -V Zi G T (x, z)V^ T (z) - T^Grfa z)^ T (z) 

OCJ e 

-T- x G t {x,z')tPt{z') 

-(Vt<h{z) + ti) [ G T {x,y)p u {y)V Zi G T {y,z) dy 



-T-^thlz) J d G T (x,y)p u (y)(G T (y,z)+G T (y,z')) dy 
+Gt(x, z)4> t {z) + Gt{x, z')(j>T{z'), 



and 
(3.35) 

< 

where 
(3.36) 



sup 



dipT 



dujf 



(0) 



gr(*)(\ViMz)\ + T- l \-ip T {z)\ + MT)(i + \M*)\ + IM-OD) 

+(1+\M*)\ + \M*')\) I gT(y)(\V Zl G T (y,z)\+T- 1 g T (y-z)) dy, 



u d {T) = 



T for 

VT for 

lnT for 

1 for 



d=2, 
d = 3, 
d = 4, 
d > 4. 



MONTE-CARLO APPROXIMATION OF HOMOGENIZED COEFFICIENTS 17 

Starting point is again the Green representation formula 

ip T (x) = / G T (x,y)p a] (y)(j)T(y) dy, 

associated with p,14j) in the form 

T~ 1 ptp T - V* • AVip T = pc/> T . 

Differentiated w. r. t. u e it turns into 

dip T {x) f dG T (x,y) ,,,,,, 
— ^ = / ~ PUvmiy) dy 

+ [ G T ( X ,y)^My)dy + [ G T ( X ,y) P Uy)^ dy. 

Jjd 0U! e Jjd OU! e 

Combined with (|3.21[) . (I3.26[) . and the Green representation formula itself, this 
shows ([3~34|) . 

We now turn to (|3.35l) and treat each term of the r. h. s. of (|3.34[) separately. 
We begin with the third line of (|3.34[) . appeal to (I3.32[) . then bound the gradient 
of the approximate corrector (Vi^T^)! by the approximate corrector |</>t(z)I + 
|</f>T(z')l itself, control the Green's function Gt by gr, and use (I3.22[) to estimate 
the supremum in uj e of the gradient of the Green's function |V Zi Gr(y, x)\. This 
term is thus controlled by the second term of the r. h. s. of (|3.35[) . The term in the 
fourth line of (|3.34l) is also estimated by the second term of the r. h. s. of (|3.35[) 
using (|3.27D (and the uniform bounds |VGt|, T~ 1 / 2 Gt < 1), whereas the last two 
terms of (|3.34[) are bounded by the hrst term of the r. h. s. of (|3.35[) using (|3 . 2T[) . 
The subtle terms are the first three ones, for which we have to estimate the suprema 
of \Viip T (z), \iPt{z)\, and \ip T (z')\ w. r. t. u e . 

We begin with the following two estimates 

(3.37) su P |Vt(*)| < \Mz)\ + (\Mz)\+\Mz')\ + l)MT) 

+ sup \Viip T {z)\, 

(3.38) supIV^tWI < \ViMz)\ + {\Mz)\ + \Mz')\+l)vd(T) 

+T- 1 sup|^ T (z)|, 

which — considered as a system of two coupled ODEs — show that there exists 
some T* > such that for all T ^ T* , 

(3.39) su P |Vt(z)| < |Mz)l + + l<Mz')l + iWCO 

+|V<Vt(*)|, 

(3.40) suplViVr^)] < \V t Mz)\ + (\Mz)\ + \Mz')\ + l)MT) 

+T- l \^ T {z)\. 

To prove (|3.37p we consider p.34|) as an ODE on iPt{z), bound i/jt(z') by iPt(z) + 
\Viil>T(z)\, and use that the four last terms of the r. h. s. of p.34[) are bounded by 
the second term of the r. h. s. of (|3.35p . as discussed above. Hence (|3.34p turns 
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into 



duj t 



< sup{|V i4 G r (a;, z)\\\7^ T (z)} + T^G^z, z)\M*)\ 
+T- 1 G T (^, z')(\tPt(z)\ + sup |Vi^r(«)) 

+{i + \Mz)\ + \Mz')\) 

x / 3 T(y)(|V 2! G T (y^)|+T-V(2/-2)) dy. 



Using the uniform bounds \VGt\,T~ 1 Gt % 1, and replacing the gradient of the 
Green's function by the Green's function itself in the integral — which we then 
control by Vd{T) — , this estimate yields (|3.37l) by integrating the ODE. We now 
turn to (|3.38[) and infer from (|3.34l) that 



-W Zi W Zt G T (z, z)Viih{z) - T~ l V Zi G T {z, z)Mz) 

-(V^r(js)+&) / V z ,G T (z,i/)fc(i;)V 2 ,GT(i/,z) dy 

-rVW / V Zi G T (z,y)p u! (y)(G T (y,z) + G T (y,z')) dy 
Jz d 

+V 2i G T (z, z)<j> T (z) + V Zl G T (z, z')<j> T {z'). 



Proceeding as from (|3.34p to (|3.37p . this implies (|3.38p . and therefore (|3.39|) and 
lpT30|) . Combining the inequality < \^t{z)\ + \V l tp T {z)\ with ([3~^)l and 

(|3.40|) yields the last estimate we need: 



(3.41) sup \Mz')\ S \M*)\ + (\Mz)\ + \<h{z')\ + l )MT) + |ViVr(*)|. 



We are finally in position to conclude the proof of (|3.35p : the first three terms 
of the r. h. s. of (pQ4"]) are estimated by (HOP)) . f3T55| . and (j3"Hj) . Replacing then 
the gradient of the Green's function by the Green's function itself, and the Green's 
function by gx, we obtain (|3.35p . as desired. 

Step 5. Proof of (|3~T8| . 

This is an application of the variance estimate [GOlOal Lemma 2.3] on ipx based 
on (|3.35p . In particular, 



(3.42) 




We distinguish the contributions of the two terms of the r. h. s. of (|3.35p in this 
sum and use the notation 

A e := 3t (z)(|V^tW|+T-VtW| + ^(T)(1 + |0tW| + |0t(^')I)), 

B e := (l + |0 T (z)| + |0 T (z')|) / .9T(y)(|V 2 ,G T (2/,z)|+r- 1 .9T(y- Z )) dy. 
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The contribution of the first term is estimated as follows: 



■ feE 



< 



]T <. 9T (z) 2 (|v^ T (z)| 2 + t- 2 \m*)\ 2 + MT)\i + \M*)\ 2 + \Mz')\ 2 ))) 



I'ti 



I Y.9T{zf (<|V^| 2 )+T- 2 (V4)+^(T) 2 (l + <4))) 
zei d J 

< v d {T) ( (| VV^t| 2 )+T- 2 (^)+^(T) 2 (1 + (#))), 

by stationarity of tpT, ipT, and VV>t- Combined with (|3.33[) for q — 2, and (|3.17p 
(which is proved below), this turns into 



(3.43) 



< 



T 3 In 9 T 


for 


d = 2, 


T 3/2 


for 


d = 3, 


T\nT 


for 


d = 4, 


T 


for 


d > 4. 



Estimate (|3.17[) is a consequence of the estimate (|3.16D of the spectral exponents 
and of the spectral representations 

1 



m = 



de D (A) 



«<|v^t| 2 ) 



(T- 1 +A) 4 
Jo A2 de o( A ) + J J de (A) 

de D (A) 



A 



(T- 1 +A) 4 



r 1 i r 00 i 

T jf — d eo (A)+ j - de B (A), 



since J °° j de B (A) < 1, and by definition (|2.6[) of the spectral exponents. 

We now turn to the second term of the r. h. s. of (I3.35[) . and note that it coincides 
with the term treated in [GM101 Step 3, Proof of Lemma 5] so that we have 



(3.44) 



( B e) 



< 



eel 



T 2 In 9 T 


for 


(2 = 2, 


T 3/2 


for 


d = 3, 


T 


for 


<2 = 4, 


VT 


for 


d = 5, 


InT 


for 


d = 6, 


1 


for 


d > 6. 



Hence, (f3T42]) . (j3~35|) . (j3^3| and (j3T44l) yield (j3~18| for d > 2 using (j3~12)) and 
(|3.15[) . Note that (|3.43p is not optimal for d — 2 in view of the estimate of (4>t) 
in (|3.17p (which is sharper). This comes from our estimate of sup w iPt(z) where 
we have replaced a gradient of the Green's function by the Green's function itself 
and therefore given up some potential better decay (at least integrated on dyadic 
annuli). In high dimensions however, there is no loss because <?t is then square- 
integrable itself. Anyway, the optimal spectral exponents for d = 2 have been 
already obtained in (|3.16p . 

Step 6. Proof of (13^41) . 

We quickly show how (|3.18p allows to get the sharper spectral exponents p. 41) . 
This is the same argument as in Step 5, except that (|3.17p can now be boostrapped 
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to using ([XTg]) and 



t 



(VtpT ■ AVip T ) = 



(T-i + A)* d6D(A) 



^3 de (A)+Ti i 



(T-i + A)* d6D(A) 



— de (A)+T 2 



- de (A) 



t de D (A), 



where T* < oo has been defined in Step 4. Proceeding as in the end of Step 5 this 
yields the spectral exponents (I3.4[) . and concludes the proof of the proposition. □ 

4. The random fluctuations 

In this section, we show that the computable quantity A n (t) defined in (jl.5l) is 
a good approximation of a 2 , in the sense that its random fluctuations are small as 
soon as n/t 2 is large. We write N* for N \ {0}. 

Theorem 4.1. There exists c > such that, for any n G N*, e > and t large 
enough, 

' 2 

ne 



\A n {t) -a 2 \ > e/t 



^ exp 



In order to prove this result, we rewrite A n (t) as A n {t)/p n , where 
p(uM)(£ ■ Y^{t)) 2 + ■■■+ P (lo^){£ ■ Y^(t)) 2 



A n (t) = 



and 



Pn 



ntE\p] 
nE[p] ' 



Both are sums of independent and identically distributed random variables un- 
der P®, thus enabling us to use standard tools from large deviations theory. We 
start with the following classical but instructing fact. 

Proposition 4.2. There exists cq > such that, for any n £ N and any small 
enough e > 0: 



\Pr, 



1 ^ e] < exp (— ne 2 /co) 



Proof. Let p(oj) = p(u>) — E[p]. It suffices to show that, for some C > 0, 



(4.1) 



\p(uj {1) ) + ---+p(Lu (n) )\ ^ne] sC exp(-ne 2 /C)- 



Chebyshev's inequality implies that, for any A 0, 



p(LO^) + ---+p(J^) > 



< e 



^ e 



— nsX-n 



exp(A(p(w«) + ---+J^ w )) 



-ne\ 



E [exp(ApM)] n . 



Using a series expansion of the exponential, one can check that there exists C > 
such that, for any A small enough, 

(4.2) lnE[exp(Ap(w))] < CA 2 . 

As a consequence, for any A small enough, 

P® fp(w (1) ) + • • • + p^W) > el < exp (-n(Ae - C*A 2 )) , 
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and for A = ej (2C), the latter term becomes exp(— ne 2 /(4C)). The event + 
■ ■ ■ + p(u/™)) ^ — ne can be handled the same way, and we thus obtain (|4.ip . □ 



What makes the proof of Proposition 14.21 work is the observation (|4.2p that the 
log-Laplace transform of p — p — E[p] is quadratic close to the origin. In order to 
prove the corresponding result for A n (t), we will need to control the log-Laplace 

transform of ^' Y ^ a\ uniformly in t. To that end, we use a sharp upper bound 

on the transition probabilities of the random walk recalled in the following theorem. 
We refer the reader to |HS93j or jWol Theorem 14.12] for a proof. 

Theorem 4.3. There exists a constant c\ > such that, for any environement w 
with conductances in [a,f3], any icN* and x G Z d , 

From Theorem 14.31 we deduce the following result. 
Corollary 4.4. Let c\ be given by Theorem \4.3\ For all A < 1/ci, one has 

\Y(t)\ 2 ' 



sup Eo 

teN' 



exp A 



< +00. 



Proof. Let S = 1/ci - A. By Theorem Ol 



En 



,X\Y(t)\*/t 



If the sum ranges over all x € (N*) d , it is easy to bound it by a convergent integral: 
t -d/2 J- e-^l 2 /' < t~ d ' 2 f e- 5 l*l 2 /' dx = f e~ s ^ 2 dx. 

By symmetry, the estimate carries over to the sum over all x € (1*) d . The same ar- 
gument applies for the sum over all x — {x\, . . . , Xd) having exactly one component 
equal to 0, and so on. □ 



The following lemma contains the required uniform control on the log-Laplace 

t 



transform of } ll> o\ 



Lemma 4.5. There exist Ai > and C2 such that, for any A < Ai and any teN* 



hiE 



exp A 



t 



<c 2 A 2 



Proof. It is sufficient to prove that there exists C3 such that, for any A small enough 
and any t, 



En 



exp A 



t 



< 1 + c 3 a 2 



We use the series expansion of the exponential to rewrite this expectation as 



+°° \fc 



fc=0 



(i-Y{t)) 2 J* 
1 Gt 



The term corresponding to k — is equal to 1, whereas the term for k = 1 van- 
ishes. The remaining sum, for k ranging from 2 to infinity, can be controlled using 
Corollary 14.41 combined with the bound 



E 



^ 2E 



(±Y(t)) 
t k 



2k 



which follows from the definition of of and Jensen's inequality. 



□ 
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We are now in position to prove Theorem 14. II 
Proof of Theorem \4-l\ Starting point is the inequality 

(4.3) Pf [A n (t)-oj >2e/t 

< Pf [A n (t) - a\ > e/t) + Pf [A^ip- 1 - 1) > e/t] . 

We treat both terms of the r. h. s. separately. For the first one, the key observation 
is that (recalling the definition of P given in (jl.lt ) 

'(£-yW(i)) 2 + --- + (£-^ (n) W) 2 



P® [A n {t) - o\ > e/t] = P® 



nt 



o\ > e/t 



Let A > 0. As in the proof of Proposition l4~2l we bound this term using Chebyshev's 
inequality: 

(4.4) 

P®[A n (t)-a?>e/t] 

•(£-r«(t)) 2 + ... + (£ •yW(t)) 2 



«J En 



exp A 



exp A 



exp 



nXe 



exp 



nXe 



By Lemma T4. 5 1 the r. h. s. of (14.4[) is bounded by 

Ae 
t 



exp \n I C2A — 



for all A small enough. Choosing A = e/2c2t (which is small enough for t large 
enough), we obtain 

(4.5) Pf [A n (t) - <j\ > e/t] < exp , 

as needed. 

We now turn to the second term of the r. h. s. of (|4.3|) . From inequality (|4.5p . 
we infer that there exists M > such that 

-2 



t [A n {t) > M] < exp - 



ne 
4^ 



Since p n is almost surely bounded by a constant, it is enough to evaluate the 
probability 

Pf\P?-l>e/(Mt)], 
which is controlled by Proposition ^. 21 

We have thus obtained the required control of the 1. h. s. of (|4.3[) . The probability 
of the symmetric event 

P® [a 2 - A n (t) > 2e/t 
can be handled the same way. □ 

5. Numerical validation 

In this section, we illustrate on a simple two-dimensional example the sharpness 
of the estimates of the systematic error and of the random fluctuations obtained in 
Theorems O and O 

In the numerical tests, each conductivity of B takes the value a = 1 or (3 = 4 with 
probability 1/2. In this simple case, the homogenized matrix is given by Dykhne's 
formula, namely Ab om = y/aflld = 2Id (see for instance [G1101 Appendix A]). For 
the simulation of the random walk, wc generate — and store — the environment 
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t 


10 


20 


40 


80 


160 


320 


640 


K(t) 


10 5 


3000 


3000 


1000 


500 


100 


20 


Systematic error 


1.27E-01 


7.43E-02 


4.17E-02 


2.46E-02 


1.26E-02 


6.96E-03 


3.72E-03 



Table 1. Systematic error in function of the final time t for K(t)t 2 realizations. 




Figure 1. Systematic error in function of the final time t for 
n(t) = K(t)t 2 realizations 

along the trajectory of the walk. In particular, this requires to store up to a 
constant times t data. In terms of computational cost, the expansive part of the 
computations is the generation of the randomness. In particular, to compute one 
realization of A t 2 (t) costs approximately the generation of t 2 x At = 4£ 3 random 
variables. A natural advantage of the method is its full scalability: the t 2 random 
walks used to calculate a realization of A t i{t) are completely independent. 

We first test the estimate of the systematic error: up to a logarithmic correction, 
the convergence is proved to be linear in time. In view of Theorem 14.11 typical 
fluctuations of t(A n ^(t) ~ a 2 ) are of order no greater than t/ yjn{t), and thus 
become negligible when compared with the systematic error as soon as the number 
n(t) of realizations satisfies n(t) ^ t 2 . We display in Table [T] an estimate of the 
systematic error obtained with n(t) — K(t)t 2 realizations. The systematic error 
is plotted on Figure [T] in function of the time in logarithmic scale. The apparent 
convergence rate (linear fitting) is —.85, which is consistent with Theorem 13.11 
which predicts —1 and a logarithmic correction. 

We now turn to the random fluctuations of A n ^ (t). Theorem 14.11 gives us a large 
deviation estimate which essentially says that the fluctuations of t(A n ^(t) — a 2 ) 
have a Gaussian tail, measured in units of tj \/n(t). The Figures [H[5] display the 
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500 - 



550 - 



500 - 



150 - 



00 - 



50 - 








Figure 2. Histogram of the rescaled fluctuations for t = 10 



350 




Figure 3. Histogram of the rescaled fluctuations for t = 20 



histograms of i(A t2 (f) - of) for t = 10,20,40 and 80 (with 10000 realizations of 
A t 2(t) in each case). As expected, they look Gaussian. 
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Figure 4. Histogram of the rescaled fluctuations for t = 40 
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Figure 5. Histogram of the rescaled fluctuations for t = 80 
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